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Foreword 

The  digital  implementation  of  the  segmented  tire  model  is  an  exten¬ 
sion  of  the  study  reiorted  in  the  basic  report  and  was  developed  at  the 
U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  as  a  part  of  the 
vehicle  mobility  research  program  under  DA  Project  1T062103A046,  "Traffic- 
ability  and  Mobility  Research,"  Task  03,  "Mobility  Fundamentals  and  Modei 
Studies,"  sponsored  by  the  Directorate  of  Development  and  Engineering, 

U.  S.  Army  Materiel  Command. 

The  study  described  herein  was  accomplished  by  personnel  of  the  Ve¬ 
hicle  Dynamics  Section  (VDS),  Mobility  Research  Branch,  Mobility  and  En¬ 
vironmental  Division,  under  the  general  supervision  of  Messrs.  W.  G. 
Shockley  and  S.  J.  Knight,  and  under  the  direct  supervision  of  Dr.  D.  R. 
Freitag  and  Mr.  A.  J.  Green.  Mr.  N.  R.  Murphy,  Jr.,  of  the  VDS  and 
Mr.  J.  F.  Smith  of  the  Mathematics  and  Analysis  Section,  Electronic  Com¬ 
puter  Branch,  wrote  the  computer  programs.  Mr.  Murphy  prepared  this 
appendix. 

COL  Levi  A.  Brown,  CE ,  was  Director  of  WES  during  the  preparation  of 
this  appendix.  Mr.  F.  R.  Brown  was  Technical  Director. 
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Conversion  Factors,  British  to  Metric  Units  of  Measurement 


British  units  of  measurement  used  in  this  report  can  be  converted  to  metric 


units  as  follows: 

Multiply 

By 

To  Obtain 

inches 

2.5U 

centimeters 

pounds 

4.444 

newtons 

pounds  per  square  inch 

0.689476 

newtons  per  square  centimeter 

IX 


Summary 


This  appendix  presents  the  procedures  for  digital  implementation  of 
the  segmented  tire  model,  developed  in  the  basic  report  for  an  analog  com¬ 
puter.  Two  procedures  are  required:  (a)  determination  of  the  segment 
spring  coefficient  from  measured  load-deflection  test  results,  and  (b)  com¬ 
putation  of  the  resultant  force  vector  transmitted  to  the  axle.  Digital 
programs  for  both  procedures  were  written  in  Fortran  IV  for  a  GE-U20  sys¬ 
tem,  and  are  included. 
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DYNAMICS  OF  WHEELED  VEHICLES 


A  MATHEMATICAL  MODEL  FOR  THE  TRAVERSAL  OF  RIGID 
OBSTACLES  BY  A  PNEUMATIC  TIRE 

APPENDIX  B:  DIGITAL  IMPLEMENTATION  OF 
SEGMENTED  TIRE  MODEL 

1.  For  digital  implementation  of  the  segmented  tire  model,  basi¬ 
cally  two  procedures  are  required:  (a)  determination  of  the  segment  spring 
coefficient  from  measured  load-deflection  test  results,  and  (b)  computation 
of  the  resultant  force  vector  transmitted  to  the  axle.  Therefore,  two 
computer  programs  have  been  prepared. 

Program  1 

2.  The  first  program  (fig.  Bl)  is  used  to  compute  the  average  de¬ 
flection  of  each  segment  for  a  given  vertical  center-line  deflection.*  It 
is  written  in  Fortran  IV  and  programmed  to  run  on  a  GE-420  system.  The 
only  input  requirement  is  a  card  containing  the  tire  radius  in  inches,  the 
number  of  segments  for  180  deg  (current  status  of  the  program  limits  the 
maximum  number  of  segments  to  6o),  and  the  desired  maximum  vertical  center- 
line  deflection  in  inches  multiplied  by  10.  This  maximum  deflection  serves 
as  a  cutoff  criterion  for  the  computer  and  should  be  a  number  that  does  not 
exceed  the  tire  section  height.  The  actual  maximum  deflection  is  multi¬ 
plied  by  10  to  ensure  an  integer  value.  This  is  a  restriction  imposed  by 
the  "DO-loop"  that  successively  increments  the  center-line  deflection. 

3.  A  representative  computer  printout  is  shown  in  fig.  B2  to  illus¬ 
trate  the  format  of  the  output  and  to  demonstrate,  in  conjunction  with  the 
following  example,  how  the  results  are  used  to  determine  an  appropriate 
segment  spring  coefficient: 


*  Appendix  A  to  this  report  presents  procedures  for  computing  the  neces¬ 
sary  effective  radial  deflections. 
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a.  Problem.  A  9*bO“l4,  2-PR  tire  with  a  l4-in.*  radius  i  s  in¬ 
flated  to  30-psi  prer.r.ure;  at  this  pressure  it  deflects 
.1  In.  under  a  yuO-.lb  load.  Divide  the  lower  half  ( l8'D  leg) 
of  the  tire  into  l8  equal  segments  and  determine  the  segment 
spring  coefficient. 

h.  Solution.  Prepare  an  input  data  card  containing  the  tj re 
radius  in  inches  (.14),  the  number  of  segments  in  180  deg 
( 18),  and  the  maximum  vertical  center-line  deflection  in 
in-hoc  multiplied  b,y  10  (say  3  x  10  =  30)-  Run,  using  the 
program  »n  fix.  Rl.  Since  the  coordinate  (i.O,  YbO)  was 
■hosen  from  the  load-deflection  curve,  examine  the  computer 
printout  (fig.  B2)  and  find  the  values  of  segment  deflections 
for  a  center-line  deflection  of  1  in.  Two  segments  on  each 
side  of  the  vortical  center  line  are  seen  to  be  influenced 
at  this  particular  deflection.  Compute  the  segment  spring 
coefficient  K  from  'the  enuation** 


2  5]  A.  cos  0. 
i=l  1  1 


=  _ _ 760  _ 

2(0. 9  A  cos  +  0.3SU  cos  15c) 

_  760 
2.5b 

=  297  lb/ in. 

This  value  is  assigned  as  the  spring  coefficient  for  each 
segment  and  will  change  only  if  the  number  of  segments  or 
the  inflation  pressure  is  changed.  Choosing  a  small  deflec¬ 
tion,  such  as  0.2  or  0.4  in.,  as  a  basis  for  determining 
the  spring  coefficient  would  most  likely  yield  a  K-value 
somewhat  different  from  one  with  a  larger  value  of  deflec¬ 
tion  as  a  basis.  Comparisons  of  measured  load-deflection 
curves  with  those  computed  for  K-values  obtained  from  the 
segmented  tire  model  at  various  deflections  revealed  that 
greater  accuracy  and  consistency  were  obtained  when  K  was 
computed  from  a  deflection  of  1  in.  This  deflection  is 
generally  in  the  area  of  maximum  curvature  on  the  load- 
deflection  curve. 


v  f  table  of  factors  for  converting  British  units  of  measurement  to 
meiric  unite,  is  presented  on  page  ix. 
y-  See  equation  development  in  paragraph  21  of  the  main  report. 
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4.  The  second  program  (fig.  B3)  is  used  to  compute  the  vertical  and 
horizontal  components  of  the  resultant  force  vector  transmitted  to  the 
axle.  It,  too,  is  written  in  Fortran  IV  for  the  GE-420  system.  It  is  a 
complete,  self-contained  program  that  was  assembled  from  a  section  of  a 
more  comprehensive  one  that  represented  a  multidegree-of-freedom  vehicle. 
Taken  out  of  context,  it  contains  none  of  the  dynamics  of  the  problem,  i.e. 
the  influence  of  the  coupled  differential  equations  describing  the  dynamics 
of  the  sprung  mass  of  the  system  has  been  removed.  Certain  modifications 
were  made  regarding  the  input  requirements,  so  the  program  demonstrates, 

in  a  static  sense,  how  the  components  of  the  resultant  wheel  spring  force 
are  computed  with  the  segmented  tire  concept.  With  other  slight  modifi¬ 
cations,  this  program  can  be  adapted  as  a  subprogram  to  describe  the  tire 
compliance  in  other  vehicle  dynamics  programs. 

5.  The  solution  is  based  on  a  purely  geometric  approach  that  treats 
each  discrete  axle  movement  interval  in  terms  of  space-oriented  coordinates 
of  the  terrain  and  wheel  center  with  respect  to  a  fixed  reference  frame. 

The  current  program  will  handle  up  to  100  terrain  profile  points  and  up  to 
24  tire  segments  equally  divided  about  a  vertical  axis  through  the  axle  of 
each  wheel.  Experience  has  indicated  that  twelve  10-deg  segments,  six  on 
each  side  of  the  vertical,  are  sufficient  to  describe  the  tire  compliance 
for  most  terrain  conditions  and  obstacle  configurations.  The  angle  to  each 
segment  center  line  is  measured  from  the  vertical  and  is  considered  posi¬ 
tive  in  a  counterclockwise  direction. 


Computer  calculations 


u.  The  following  sequence  of  calculations  is  used  in  this  second 


program : 

a.  Dimension  the  appropriate  space  for  segment  angles,  segment 
spring  coefficients,  terrain  profile  points,  segment 
forces,  etc. 

_b.  Provide  a  "degree-to- radian"  converter  =  Jt/l80. 

c.  Read  in:  Coordinates  of  wheel  center  (Xl,Yl),  tire  radius, 
number  of  tire  segments,  number  of  terrain  points,  number  of 
positions  to  be  calculated,  and  X  and  Y  increments. 
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Check  for  end.  of  data;  if  no  more  data,  exit  computer. 

Read  in  and  store  all  segment  angles,  spring  coefficients, 
and  terrain  points. 

Print  headings  and  terrain  profile  table. 

Designate  a  variable  name  for  the  number  of  tire  segments; 
e.g.  NOT  =  NSEG. 

Determine  maximum  horizontal  projection  of  tire,  i.e. 

XRET  =  XI  -  WRAD,  XFET  -  XI  +  WRAD. 

Initialize  ITR  =  2. 

Change  all  segment  angles  from  degrees  to  radians. 

Compute  the  coordinates  of  each  (undeflected)  tire  segment 

(1)  location  with  respect  to  the  wheel  axle: 

CX(i)  -  SIN(  DSEG(  i )  *WRAD) 

CY(  i  )  =  -COS(DSEG(  i)  *WRAD) 

(These  then  remain  constant  for  a  given  wheel.) 

Set  VW(1)  =  0  VW(2)  =  0. 

Now  determine  the  coordinates  of  each  (undeflected)  segment 
with  respect  to  the  fixed  reference  system,  e.g. 

KXY(l,i)  -  XI  +  CX( i ) 

RXY(2,i)  =  Y1  +  CY(i) 

At  this  time  compute  intersections  of  the  tire  segment 
centers  with  terrain.  For  each  tire  segment: 

(  l)  Set  DMIN  -  J.E10,  i.e.  some  large  number  definitely 
greater  than  tire  radius. 

(2)  Compute  slope  of  tire  segment  (i),  e.g.  SM  ~  CY(i)/cX(i) 

(j)  KTR  -  ITR.  This  is  a  counter  to  account  for  terrain 
stations. 

(4)  Set  X2  =  0,  i.e.  first  terrain  point  begins  at  X  =  0. 

(5)  Recall  Y2  from  proper  storage  location;  now  there  is 
terrain  coordinate  X2,Y2. 

(t)  Recall  X3,Y3,  i.e.  the  next  forward  terrain  coordinate. 

(7)  Obtain  the  slope  of  this  line  (terrain  segment). 

(8)  Compare  it  with  the  slope  of  i^*1  tire  segment. 

(9)  If  slopes  are  not  the  same,  the  lines  will  intersect; 
use  the  point  slope  method  to  obtain  coordinate  of  inter 
section,  i.e.  XX, YY  for  i^h  tire  segment.  (This  is 


coordinate  of  deflected  tire  segment  with  respect  to 
fixed  reference  system.) 

(10)  Nov/-  determine  whether  this  coordinate  falls  within  the 
tire  circumference. 

(11)  If  not,  then  increment  the  terrain  segment  and  recheck 
as  before  until  a  coordinate  is  determined  that  falls 
within  the  tire  circumference. 

(12)  Comp u be  the  radial  distance  of  this  intersection  from 
the  wheel  center  and  set  this  distance  equal  to  the 
variable,  DMIN.  This  represents  a  new  "minimum"  length. 

(13)  Check  this  length  against  last  minimum  value. 

(lit)  If  it  is  smaller,  let  the  coordinates  of  this  deflected 
tire  segment  be:  RXY(l,i)  =  XX;  RXY(2,i)  -  YY. 

(15)  If,  however,  it  is  larger,  check  to  see  whether  forward 
edge-  of  tire  is  ahead  of  the  terrain  station. 

(lo)  If  so,  then  advance  the  station  by  one.  Go  back  to 
step  ( 1 ) . 

(1?)  If  not,  go  to  next  tire  segment  and  repeat  the  proce¬ 
dures  described  above  until  all  intersections  of  tire 
segment  centers  with  terrain  segments  have  been 
computed. 

This  operation  should  locate  the  positions  (with  respect  to 
the  fixed  and  moving  reference)  of  each  deflected  tire 
segment . 

Nov/  compute  the  actual  deflection  of  each  tire  segment  be¬ 
ginning  as  follows: 

RXY( 1 , i )  =  XI  -  RXY(l,i) 

RXY(2,i)  =  Y1  -  RXY(2,i) 

These  are  vertical  and  horizontal  components  of  new  segment 
length. 

Store  these  values  in  a  temporary  storage  M(l),AA(2). 

Compute  the  radial  distance  of  each  segment  using  the 
Pythagorean  Theorem,  i.e.  TMAG( i )  =  square  root  of  sum  of 
squares  of  the  respective  components. 

Obtain  the  sin  and  cos  of  each  segment  angle  as  follows: 

RXY(  l,i)  =  AA(  1  )/TMAG( i )  =  sin  of  the  segment  angle 
RXY(  2 ,  i )  -  AA(  2  )  /TMAG(  i  )  =  cos  of  the  segment  angle 

Compute  the  radial  force  in  each  segment,  i.e. 
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TMAG(  i )  =  ( WRAD  -  TMAG(  i  )  )*  SEGK(  i )  =  K6 


where  6  =  WRAD  -  TMAG(  i ) 

t.  Compute  the  horizontal  and  vertical  components  of  liio  re¬ 
sultant  force  vector,  e .  g . 

VW(l)  -  ITMAG{  i  )  *  RXV(l,i)  -»  horizontal  component 
VW(2)  =  LTMAG(  i )  *  RXY(2,i)  -»  vertical  component 

u.  Print  results. 

v.  Increment  tiro  center  j)Osition  and  repeat  procedures. 

7.  Computation  of  resultant  tire  spring  force  vector  is  summarized 
as  follows: 

a.  The  points  where  all  tire  segment  centers  (spring  locations) 
intersect  the  undeflected  tire  circumference  are  located. 

b.  The  points  where  tire  segment  centers  intersect  the  terrain 
are  located. 

c.  The  radial  reduction  in  length  of  each  tire  segment  center 
intersecting  the  terrain  is  computed. 

d.  The  force  vector  for  each  segment  is  computed  using  the  dis¬ 
placement  of  the  segment  center  and  the  segment  spring 
coefficient. 

e.  The  segment  force  vectors  are  summed  to  yield  the  vertical 
and  horizontal  components  of  the  resultant  spring  force  vec¬ 
tor  for  the  wheel. 

Symbols  used  in  the  second  program  are  defined  below: 

CX,CY  Coordinates  (with  respect  to  the  wheel  axle)  of 
points  where  segment  center  lines  intersect  the 
undetected  tire  circumference 

DSEG  Segment  angle  in  degrees;  the  angle  measured  from 
a  vertical  line  through  the  wheel  axle  to  center 
line  of  segment  (positive  value,  counterclockwise ) 

DTR  Degree-to-radian  conversion  factor 

HGT  Elevation  of  profile 

NPOS  Number  of  wheel  positions  to  be  computed 

NSEG,NWN  Number  of  tire  segments 

NTP  Number  of  terrain  points  to  be  used 


Notation 

8. 
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RXY  Absolute  coordinates  (with  respect  to  fixed  X-Y 
reference)  of  the  points  where  segment  center 
lines  intersect  the  undeflected  tire  circumference 

SEGK  Segment  spring  coefficient  (can  vary  with  each 
segment) 

ST  Slope  of  a  terrain  segment 
STA  Horizontal  distance  of  terrain  profile 
TMAG  Radial  force  of  segment 

VW(l)  Horizontal  component  of  force 

VW(2)  Vertical  component  of  force 
WRAD  Wheel  radius 

XX, YY  Absolute  coordinates  of  the  intersection  points 
of  the  terrain  with  tire  segment 

XI, Y1  Coordinates  of  wheel  center 

X2,Y2  Absolute  coordinates  of  the  rear  terrain  station 

X3,Y3  Absolute  coordinates  of  the  adjacent  terrain  sta¬ 

tion  in  increasing  X-direction 

XFET,XRET  Forward  and  rear  extremes  of  tire  (horizontal) 

XINCR, YINCR  Increment  of  axle  movement  of  wheel  in  X-  and 
Y-directions 


Summary 

9.  The  required  input  information  for  the  second  program  consists 
of  (a)  the  coordinates  of  the  tire  center,  (b)  the  undeflected  tire  radius, 
(c)  the  number  of  segments  to  be  used,  (  d)  the  number  of  terrain  coordi¬ 
nates,  (e)  the  number  of  axle  positions  to  be  calculated,  (f)  the  X-  and 
Y-increments  (i.e.  the  movements  of  the  axle),  (g)  each  segment  angle  in 
degrees,  (h)  each  segment  spring  constant  (may  be  different  for  each  seg¬ 
ment  if  desired),  and  (i)  the  terrain  profile  coordinates.  The  tire-center 
coordinates  are  the  only  inputs  that  do  not  remain  constant.  When  the  pro¬ 
gram  is  used  in  an  appropriate  vehicle  dynamics  program  that  describes  the 
motions  of  the  sprung  masses,  these  coordinates  are  determined  at  each  time 
step  by  the  differential  equations  that  describe  the  dynamics  of  the  ve¬ 
hicle  system.  The  coordinates  (Xl,Yl)  of  the  tire  center  are  used  to  com¬ 
pute  forward  and  rear  extremes  of  the  tire  for  each  increment  of  movement. 


This  is  accomplished  by  subtracting  or  adding  the  tire  radius  (WRAD)  to 
the  coordinate  XI.  (See  fig.  b4.)  The  program  can  then  be  used  to  compute 
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the  location  of  the  points,  relative  to  the  axle,  where  each  segment  cen¬ 
ter  intersects  the  tire  circumference.  These  values  never  change  and  are 
stored  in  the  computer  under  the  identification  CX(i)  and  CY{  i ) . 

10.  Then  for  each  increment  of  horizontal  axle  movement  (XINCR), 

the  program  proceeds  as  foiLows: 

The  points  of  the  wheel  at  each  segment's  center  line 
are  located  relative  to  the  fixed  (X-Y)  reference 
frame  and  stored  as  coordinates  under  the  identifica¬ 
tion  RXY(l,i)  and  RXY(2,i).  This  gives  the  complete 
orientation  of  each  undeflected  segment  "spring."  The 
program  then  proceeds  into  a  .loop  that  calculates  the 
intersections  of  the  segment  center  lines  and  the  ter¬ 
rain  profile  segments,  which  are  constructed  from 
straight-line  connections  of  the  terrain  profile  points. 

These  intersections  are  determined  by  the  point-slope 
method,  which  compares  the  slopes  of  each  segment  cen¬ 
ter  line  ( proceed ing  in  a  counterclockwise  sequence  be¬ 
ginning  with  the  rear  segment)  with  the  slopes  of  each 
terrain  segment  within  the  forward  extreme  of  the  tire. 

A  check  is  made  in  each  instance  to  determine  if  the 
intersections  occur  within  the  periphery  of  the  unde¬ 
flected  tire.  This  indicates  whether  a  particular 
segment  is  influenced  by  the  terrain  and  the  amount 
that  each  segment  deflects  in  a  radial  direction.  The 
radial  deflection  that  is  computed  for  each  segment 
spring  is  multiplied  by  the  appropriate  spring  con¬ 
stant  to  yield  a  segment  force  vector  whose  magnitude 
and  orientation  are  known.  These  individual  vectors 
are  summed  to  yield  the  vertical  and  horizontal  com¬ 
ponents  of  the  resultant  force  vector  acting  at  the 
axle.  The  axle  is  then  advanced  to  the  next  position 
and  the  process  is  repeated. 


11.  No  provision  has  been  made  to  incorporate  segment  damping 
forces  at  this  time.  The  influence  of  tire  damping  is  currently  computed 
as  a  gross  vertical  force  from  the  relative  motion  between  the  axle  and 
a  point  directly  beneath  the  axle. 
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DIRECTION  OF  NOTION 


Schematic  diagram  depicting  space  coordinates  of 
digital  segmented  tire  program 
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